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Abstract 

The need to understand and reject backgrounds in Ge-diode detector double-beta 
decay experiments has given rise to the development of pulse shape analysis in such 
detectors to discern single-site energy deposits from multiple-site deposits. Here, 
we extend this analysis to segmented Ge detectors to study the effectiveness of 
combining segmentation with pulse shape analysis to identify the multiplicity of 
the energy deposits. 
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1 Introduction 



Zero- neutrino, double-beta decay 
(PP(0u)) studies are well motivated 
on physics grounds and several re- 
cent reviews make this case. (See for 
example Refs. [1][2][3][4].) The pro- 
cess of PP(0u) may occur in certain 
even-even nuclei, where (3 decay is 
forbidden, but only if the neutrino 
is a massive Majorana particle. Dur- 
ing this process a nucleus will change 
atomic number by two units while 
emitting 2 electrons (e.g. 76 Ge — > 76 Se 
+ 2e~). Because the rate of this lep- 
ton number violating process is pro- 
portional to the square of the effec- 
tive Majorana neutrino mass there is 
strong current interest in this decay. 
With the electrons being the only light 
particles in the final state, the sum of 
their energies is mono-energetic and 
therefore a distinct signature for the 
decay. In any /^(Oz/) experiment the 
peak, if it exists, will be superim- 
posed upon a continuum of background. 
Present experimental limits indicate 
the Pj3(0v) decay rate would be very 
low even in the larger mass double 



The signature for a f3f3(Qv) in a Ge 
detector is a localized deposit of ion- 
ization of 2.039 MeV; that is, a single- 
site event. In contrast, 7 rays of a 
MeV or more tend to multiply scat- 
ter as they interact in a solid result- 
ing in a multi-site event. Furthermore, 
there are potential cosmogenic iso- 
tope decays that can populate the 
(3(3(0is) region of interest near 2 MeV 
including the 68 Ga f3 + decay and the 
60 Co 7 cascade. All of these back- 
ground processes tend to produce multi- 
site energy deposits. This fundamen- 
tal difference in the energy deposi- 
tion process creates an opportunity 
for separating the signal from back- 
ground. 

Ge detectors have been used for sev- 
eral decades in the study of f3f3(0v) 
beginning with the initial work of Ref . 
[5] culminating with recent state-of- 
the-art results from the IGEX [6] and 
Heidelberg-Moscow [7] experiments. 
Iln a coaxial Ge detector, the shape 
of the waveform will primarily de- 
pend on the radial distribution of the 
ionization produced during an event. 
Both IGEX [8] and Heidelberg-Moscow 



beta-decay experiments now being planned^] used pulse shape analysis (PSA) 



making the identification and elimi- 
nation of background essential. 

For the case of 76 Ge, the total energy 
that the 2 electrons possess is 2.039 
MeV. This energy is above most, but 
not all, the Q- values of naturally oc- 
curring radioactive isotopes. The range 
of the (3(3(0u) electrons in a Ge crys- 
tal is no more than a couple mm. 
Therefore, solid state Ge detectors 
make a high-efficiency experimental 
apparatus with excellent energy res- 
olution for the study of this decay. 



for the latter fractions of their data 
sets in order to improve their sig- 
nal to background ratio. Both exper- 
iments used large coaxial Ge detec- 
tors (~2 kg or 80 mm diameter) re- 
sulting in a significant background 
rejection capability. 

In addition to PSA, segmented Ge 
detectors will also have a capability 
to identify multiple-site energy de- 
posits. Although this technology has 
not yet been brought to bear for /3/3(Oz/) 
it has been used [10] [11] [12] to ana- 
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lyze the multiple interactions arising 
from a 7 ray impinging on a detector 
to "track" the 7-ray path. These 7 
tracking arrays require very detailed 
electric field models for each detec- 
tor in the array [13]. These models 
are then used to simulate a library 
of pulses for comparison to the mea- 
sured pulse. This tracking capability 
utilizes image-charge formation in ad- 
jacent segments to improve event re- 
construction and probably provides 
the highest resolution localization of 
events in Ge of any technique in use 
at present. 

Recently the Majorana collaboration 
[14] has proposed to use simple PSA 
in conjunction with segmentation to 
further reduce background levels in 
(3/3 (Ov) without the need for very de- 
tailed detector field models. The de- 
tector schemes considered for this pro- 
posal include segmenting coaxial Ge 
detectors along either the axial or az- 
imuthal direction. The electronic sig- 
nals from the various segments would 
provide information regarding the ax- 
ial and azimuthal distributions of ion- 
ization and therefore should improve 
on the radial-only information pro- 
vided by PSA. Although Ge detec- 
tors have been produced with a large 
number of segments (the Gretina de- 
tector has 36 segments), the maxi- 
mum number of segments isn't nec- 
essarily the optimum for (3(3(Ois). The 
cost is greater for finer segmentation 
and the additional detector contacts 
with their associated electronic com- 
ponents and cables may add to the 
background. 

The purpose of this work is to ex- 
perimentally measure the effective- 



ness of the combination of PSA and 
segmentation in reducing the back- 
ground for (3(3(0v). Although we also 
discuss the reduction of background 
via a granularity cut (i.e., requiring 
only one detector in an array observe 
an energy deposit), we focus on the 
impact of segmentation. Both seg- 
mentation and granularity contribute 
to an array's self-shielding and hence 
the identification of single-site energy 
deposits. The results of this study 
will aid in the segmentation design 
to optimally reduce background. 



2 Experimental Setup 

The data for this study were taken 
using a CLOVER [15] detector. The 
CLOVER is a close-packed array of 
four germanium detectors, each 50 
mm in diameter and 80 mm long. 
Each detector is electrically segmented 
into two parts as shown in Fig. 1. 
The CLOVER has four high-resolution, 
cold-FET energy readouts (one for 
each crystal), and three low-resolution, 
warm-FET position readouts corre- 
sponding to the left two, middle four 
and right two segments. Coincidences 
between the energy and position read- 
outs indicate which segment (s) recorded 
energy depositions. 

The major focus of PSA is to distin- 
guish multi-site from single-site en- 
ergy deposits, and therefore a pop- 
ulation of each is required for our 
study. To this end, we follow the ex- 
ample of Refs. [8] and [9] and use the 
7-ray line from 228 Ac at 1588 keV 
and the double escape peak at 1592 
keV derived from the 2614 keV 7 ray 
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Fig. 1. A schematic sketch of the layout 
of the Ge detectors in the CLOVER. 
This figure is a rendition from the 
owner's manual. 

from 208 T1. (For reference, the mean 
free path for a 1.6 (2.6) MeV 7 ray 
in solid Ge is about 4 (5) cm.) Since 
both of these isotopes are in the de- 
cay chain for 232 Th, a single source 
can be used to obtain both peaks. 
The 1.6-MeV 7 ray from 228 Ac typ- 
ically interacts about 3 times while 
depositing its entire energy in a crys- 
tal. Hence the full-energy 7-ray peak 
from 228 Ac consists almost completely 
of multi-site events. The double es- 
cape peak (DEP), conveniently, is com- 
pletely single-site because it requires 
the annihilation 7 rays from the electron- 
positron pair production escape the 
crystal. A Th source spectrum in the 
energy region near these peaks is shown 
in Fig. 2. Although both of these pro- 
cesses are required to define the PSA 
analysis, a nearly pure sample of DEP 
events can be obtained by requiring a 
coincident 1592-keV energy deposit 
in one crystal while simultaneously 
detecting one of the 511-keV 7 rays 
in a different crystal. 

All seven channels from the clover 
are read out using a pair of X Ray In- 
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Fig. 2. The energy spectrum of events 
from a Th source near 1.59 MeV. 

strumentation Associates (XIA) [16], 
Digital Gamma Finder Four Channel 
(DGF4C) CAMAC modules. These 
are 14-bit digitizers with a 40 MHz 
sampling rate. The CAMAC crate is 
connected to the PCI bus of a Dell 
Optiplex computer running Windows 
2000. The system was controlled us- 
ing the standard software supplied 
by XIA. This data acquisition soft- 
ware runs in the IGOR Pro environ- 
ment [17] and produces binary data 
files that are read in and analyzed 
using the ROOT framework. 



3 Pulse Shape Analysis 



Figure 3 shows a current waveform 
from an event at 1588 keV and there- 
fore most likely a multiple-site en- 
ergy deposit 7-ray event. Figure 4 
shows a waveform from an event at 
1592 keV and therefore most likely 
a single-site energy deposit from a 
DEP event. The waveforms are clearly 
different. We analyze the waveforms 
using the formalism developed by 
Aalseth [18]. 

We calculated three different pulse 
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Fig. 3. The current pulse waveform for 
a pulse near 1588 keV and likely to be 



an 228 Ac 7 ray. 
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Fig. 4. The current pulse waveform for 
a pulse near 1592 keV and likely to be 
a 208 T1 DEP event. 

shape parameters (Fig. 5) for each 
pulse: the width, front-back asym- 
metry, and a normalized moment. The 
pulse width is defined by the dura- 
tion between the times where the charge 
pulse rises to the baseline plus 10% 
and 90% of the pulse height. With 
the pulse mid-point defined as that 
time half way between the 10% and 
90% points, the front-back asymme- 
try is the difference between the area 
in the front and back halves of the 
mid-point of the current pulse divided 
by the total area of the pulse. Last, 
the normalized moment parameter is 
analogous to the moment of inertia if 
we treat the current pulse as a mass 
distribution (i.e., it is a sum of the 
amplitude of the current pulse mul- 
tiplied by the square of the distance 
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Fig. 5. Cartoon illustrating idealized 
pulse shape analysis parameters. The 
monotonically raising curve is a charge 
pulse, and the shaded-to-zero curve is 
a current pulse (constructed by differ- 
entiating the charge pulse) . See text for 
a detailed discussion of the calculations 
of the pulse parametes. 

from the pulse midpoint). The nor- 
malization of the normalized moment 
parameter (J n ) is chosen to minimize 
any energy dependence. The specific 
formula is: 



12 



(1) 



where ji is the current pulse ampli- 
tude at time index %. N, N m id, and 
Nq are the time indices at the pulse 
onset, midpoint, and end respectively. 
At is the duration of time between 
clock ticks of the digitizer. 



4 Data 



Once we calculate the three PSA pa- 
rameters, we form a two-dimensional 
histogram in the asymmetry and nor- 
malized moment parameters. In our 
pulse shape analysis, we have tended 
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not to use the pulse width as an anal- 
ysis parameter because it is largely 
degenerate between the 7-ray and DEP 
peak populations. Previous studies 
have found this parameter useful [18] 
and we hypothesize that our smaller 
crystals lead to less sensitivity to this 
parameter. 

We separately fit empirical functions 
to the two-dimensional histograms of 
the remaining two parameters for DEP 
and 7-ray events in a calibration sam- 
ple of our data. The fit results in two 
probability distribution functions, one 
for each class of event (single-site or 
multi-site) that depend upon the two 
critical PSA parameters. For a given 
event, the measured PSA parameters 
are used to calculate the probabili- 
ties that it is associated with a given 
class. We assigned the event to the 
single-site or multi-site class for which 
the probability is higher. This proce- 
dure results in accepting 75% of the 
DEP events but rejects 80% of the 
228 Ac 7-ray events. The nearby 7- 
ray peaks that were not used for this 
calibration have similar rejection re- 
sults. 

One caveat to the comparison be- 
tween DEP and 228 Ac events is that 
they have moderately different spa- 
tial distributions throughout the de- 
tector. The DEP events tend to be 
somewhat nearer the detector edge 
because both annihilation 7 rays must 
exit the detector. The mean free path 
for a 511-keV 7 ray in Ge is about 2 
cm, which is smaller than the crystal 
dimensions. Hence the DEP events 
tend to fall nearer the detector edge 
where the 511-keV photon has a higher 
probability to escape. In contrast the 



distribution of the 228 Ac events is de- 
fined by the 4-cm mean free path 
of the 7 ray to its initial interaction 
point. In contrast, to DEP events, 
PP(0u) events would be uniformly dis- 
tributed throughout the Ge. 



5 Simulation 

We simulated the distribution of 228 Ac 
7-ray and 208 T1 DEP events within a 
Ge crystal of the same dimensions as 
those composing the CLOVER using 
GEANT3-GCALOR [19]. We gener- 
ated r, (j) and z coordinates for the 
energy deposits of these events. Fig- 
ures 6 and 7 show the radial and 
axial distributions for the simulated 
DEP events for a source positioned 
near the front face of the CLOVER. 

Pulse Shape Analysis effectively elim- 
inates all multi-hit events in which 
energy is deposited outside of some 
"resolution" for identifying radial ion- 
ization positions. This effective reso- 
lution also depends on the multiplic- 
ity of the Compton events, improv- 
ing as the multiplicity increases (and 
therefore ensuring the radial ioniza- 
tion distribution exceeds that expected 
for single-site events). To determine 
the effective resolution, we establish 
the resolution which corresponds to 
rejecting the appropriate number of 
multi-site events using a Monte-Carlo 
distribution for the radial separation 
between the Compton sites. This de- 
pends on the type of event (photo- 
peak vs. Compton continuum) both 
in the predicted spatial distribution 
for events and in the resolution. For 
the simulated 228 Ac events, we com- 
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pared our measured survival proba- 
bility to the predicted survival clS du 
function of effective separation. From 
these results (Fig. 8) we estimate the 
effective r separation to be 3 mm. 
That is, an event with multiple en- 
ergy deposits at radii that differ by 
more than 3 mm will be interpreted 
as a multiple-site event. Figures such 
as 8 and 9 can therefore be used to 
estimate the survival probability for 
other segmentation configurations and 
PSA spatial resolution capabilities. 

The measurements indicate that about 
3% of the DEP events were rejected 
due to the segmentation cut. We de- 
termined the number of events be- 
fore and after the cut by fitting to 
the spectrum. Even so, many of these 
removed events might originate from 
7-ray interactions that form the con- 
tinuum beneath the DEP peak. If we 
assume all the removed events are 
actually DEP events that are excluded 
because they fall on the border be- 
tween the segments, however, we can 
derive an upper estimate for the size 
of the border. The border is a diam- 
eter across the detector face. There- 
fore, assuming the DEP events uni- 
formly illuminated the front face of 
the detector, a border width of ap- 
proximately 1.2 mm would encom- 
pass 3% of the DEP events. 

To understand the border region bet- 
ter, we scanned the front face of the 
CLOVER with a collimated 137 Cs 
source. The collimator was formed 
by a 2-mm hole through a 4-inch Pb 
block. We moved this source across 
the border in steps of 2 mm near the 
border, and larger steps away from 
it. We scanned across detector num- 



ber 1 moving from the left position 
region through the border region into 
the middle position region. The scan- 
ning was performed slightly off the 
diameter to avoid the central core of 
the detector. The collimator illumi- 
nated the detector with a spot size 
of 3.5 mm diameter. Fig. 10 shows a 
plot of the average fraction of the en- 
ergy observed in crystal 1 that is as- 
signed to the left or middle positions 
for events that populate both seg- 
ments. Specifically we have plotted 
(En-EMiEn+Et) where E m (£,) is 
the energy in the middle (left) posi- 
tion. We then fit a step function con- 
volved with a Gaussian to these data 
and took the width of that Gaussian 
to be the total width corresponding 
to the finite spot size plus the width 
of the segmentation border. Since the 
diameter of the spot was known, we 
could simply subtract it off to obtain 
the width of the border region. We 
calculated this border region width 
to be 1.9 mm. That is, an event within 
about 1 mm of the demarcation line 
will populate both segments. Since 
this study was done with 7 rays, most 
of the two-segment events are most 
likely true multi-site events. There- 
fore this border width estimate is also 
only an upper limit. 



6 Analysis 

Figure 11 shows the Th source spec- 
trum near the 76 Ge double-beta de- 
cay endpoint. Presumably most of the 
events in this region were from the 
Compton tail of the 208 T1 2.6-MeV 

7 ray. Such events are a good model 
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Fig. 6. The radial distribution of simu- 
lated DEP events for a source position 
in front of the CLOVER. Ro is the ra- 
dius of the detector (25 mm). 
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Fig. 7. The axial distribution of simu- 
lated DEP events for a source position 
in front of the CLOVER. h is the to- 
tal height of the detector (80 mm). The 
front face of the detector is represented 
by h=0. 

of background for double-beta decay 
due to 7 rays originating outside the 
Ge detector. Events within a full en- 
ergy peak (such as those from 228 Ac 
at 1.58 MeV) are predominantly multi- 
site as they are comprised of several 
7-ray interactions. The fraction of events 
in the Compton tail that are single- 
site, however, is appreciable. Not sur- 
prisingly, the rejection of events in 
this 2.0 - 2.1 MeV window is less 
than that for the full energy 7-ray 
peak. Here the survival probability 
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Fig. 8. The fraction of 7-ray events re- 
moved function of the effective 
radial separation as predicted by sim- 
ulation. Our measured PSA rejection 
probability indicates an effective radial 
separation of ~ 3 mm for 228 Ac photo- 
peak events with an average multiplic- 
ity of about 3. 
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Fig. 9. The fraction of 7-ray events 
removed as a function of the number 
of segments as predicted by simula- 
tion. The predicted value for 2 seg- 
ments (35%) agrees with the measured 
value (34%, see Table 1) 

for the two cuts is 30%. In addition, 
however, a granularity cut is helpful 
for such Compton scatters. Requir- 
ing that only 1 of the 4 detectors 
has a signal above a 30-keV thresh- 
old reduces the rate in this window 
by a factor of 1.89 ± 0.023 ± 0.09, 
where the first uncertainty is statis- 
tical and the second is a systematic 
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Collimator Position (mm) 

Fig. 10. A measurement of the border 
width of crystal 1 of the CLOVER. The 
average fraction of energy assigned to 
the left and middle positions is plotted 
as a function of source location. The 
curve is a fit to the data. See text for 
discussion of the collimated source. The 
absolute value of the position scale is 
arbitrary. 

uncertainty estimated by the differ- 
ence in the results found from dif- 
ferent source placements. All our re- 
sults concerning the segmentation and 
PSA analyses were calculated on data 
sets after all events failing the gran- 
ularity test were removed. 

Single-escape peak events (SEP) are 
necessarily multiple-site energy de- 
posits although typically with fewer 
interaction points than a 1.59-MeV 
228 Ac 7-ray event. Thus the survival 
probabilities for this class of events 
is between the 228 Ac 7 ray and the 
DEP values. 

Table 1 summarizes the survival prob- 
abilities for the segmentation and PSA 
cuts and for the combined cuts. The 
statistical uncertainties are determined 
from the ratios of events before and 
after each cut. The systematic uncer- 
tainties are estimated by comparing 
the surviving fractions from three dif- 
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Fig. 11. The energy spectrum of events 
from a Th source near the 76 Ge dou- 
ble-beta decay endpoint of 2039 keV. 
The top curve is the raw spectrum, next 
is the spectrum after the segmentation 
cut, next is the spectrum after the PSA 
cut, and the bottom curve is after both 
segmentation and PSA cuts. 

ferent data sets, each with the source 
at a different location with respect 
to the CLOVER. The source was lo- 
cated near the center of the front face 
of the CLOVER, and near 2 sides of 
the CLOVER. 

The PSA survival results are compa- 
rable to those found by Aalseth [18] 
for the 228 Ac 1.588-MeV line (26%) 
and for the DEP line (80%). The sur- 
vival results are also similar to those 
of Ref. [20] (20% and 62% respec- 
tively) . 

The data in Fig. 2 and Table 1 show 
that the product of the measured sur- 
vival probabilities due to pulse shape 
analysis (P r = 0.20) and that due to 
the segmentation cut (P seg = 0.66) 
is greater than the measured survival 
probability of the joint cut (P r ,seg 
0.07). This indicates that the two cuts 
are effective in removing events from 
separate, but not necessarily exclu- 
sive, subsets of the entire data set. 
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Table 1 

A summary of the survival probabilities for each data cut for several processes in 
the CLOVER. The first uncertainties are statistical and the second are systematic. 

Process Energy Segmentation PSA Seg. and PSA 

228 Ac 7 ray 1.588 MeV (66± 1.2 ± 0.7)% (20± 0.5 ± 1.0)% (7± 0.3 ± 0.4)% 

208 X1 DEP 1 592 MeV (97± 2.4 ± 1.2)% (75± 2.0 ± 2.1)% (73± 2.0 ± 4.0)% 

208 T1 g Ep 2.103 MeV (63± 1.4 ± 1.8)% (45± 1.1 ± 4.5)% (20± 0.7 ± 2.0)% 

Compton Continuum 2.0-2.1 MeV (81± 1.6 ± 2.0)% (43± 0.9 ± 3.0)% (30± 0.6 ± 2.0)% 

208 T1 7 ray 2.614 MeV (70± 2.0 ± 4.1)% (17± 0.7 ± 2.2)% (8± 0.5 ± 0.4)% 



These measured values agree with the 
simulation results plotted in Fig. 12. 

By comparing the simulation to the 
measurements, we can develop insight 
as to how the survival probabilities 
change with improvements in the spa- 
tial resolutions. Figs. 12 and 13 show 
that, as one would expect, the joint- 
cut rejection improves as either the 
radial resolution improves or the seg- 
mentation increases. For the CLOVER, 
however, the spatial resolutions are 
rather modest. The segmentation is 
coarse and the PSA is limited due to 
this modest segmentation. In more 
highly segmented detectors, PSA on 
the induced pulses from multiple seg- 
ments can greatly improve the spa- 
tial resolution. Fig. 13 indicates that 
if the detector's effective resolution 
is good in any one dimension, the 
incremental gain in background re- 
jection by improving the resolution 
in an orthogonal direction is dimin- 
ished. In Fig. 13, the joint-cut re- 
jection is shown for a fixed effective 
radial resolution (3 mm) but as the 
angular resolution changes. The plot 
shows that as the angular resolution 
changes (i.e., the number of azimuthal 
segments increases from the Clover's 
current twofold azimthal segmenta- 
tion), a value is reached where the 



joint-cut survival probability is no longer 
better than the product of the indi- 
vidual cut survival probabilities. 

This situation leads to an optimiza- 
tion opportunity. There is no fun- 
damental need to increase the seg- 
mentation beyond the point where 
Pr,seg >P S eg x Pr- Once this condition 
is reached, the added spatial resolu- 
tion has reached a point of dimin- 
ishing returns with respect to back- 
ground rejection, but it might be costly 
to implement. From Fig. 13, this con- 
dition is met for an effective radial 
separation of 3 mm and ~5 segments 
for detectors of the size used here. 



7 Conclusion 

We have measured the rejection of 
multiple-site energy deposits and the 
acceptance of single-site energy de- 
posits in a segmented Ge detector by 
combining pulse-shape analysis with 
segmentation cuts. We have compared 
these measurements to simulation to 
better understand the complementar- 
ity of the cuts. Since the two cuts 
operate on orthogonal coordinates, 
they tend to remove events from dif- 
ferent subsets of the data and there- 
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Fig. 12. The joint-cut survival probabil- 
ity (P r:Se g, lower curve) and the product 
of the individual cut survival probabili- 
ties (P r xP ses ) as a function of effective 
spatial resolution as predicted by simu- 
lation for 228 Ac 7 rays. For both curves, 
the 4> survival is determined from a 
straight-forward segmentation cut. For 
an effective radial separation of 3 mm, 
the simulation agrees with the data. 
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Fig. 13. The joint-cut survival probabil- 
ity (P r!Seg , the mostly higher curve) and 
the product of the individual cut sur- 
vival probabilities (P r xP ses ) as a func- 
tion of effective azimuthal resolution 
as predicted by simulation for 228 Ac 7 
rays. For both curves, the r survival is 
determined by assigning all events with 
a radial separation greater than 3 mm 
to be multiple site. 



fore the joint action of both tends to 
be more effective than is predicted 
from the product of each individual 
cut. However, as the spatial resolu- 
tion in one of the dimensions becomes 
better, the impact of the additional 
cut is moderated and for very good 
resolution, the joint cut is not as ef- 
fective as would be predicted by the 
product of the individual cuts. 
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